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SEEDING CRITERION FOR NONEQUILIBRIUM 
MAGNETOHYDRODYNAMIC GENERATORS* 
by Frederic A. Lyman and Eli Reshotko 
Lewis Research Center 

SUMMARY 

The various consequences of the addition of seed to a carrier gas in a magnetohy- 
drodynamic (MHD) generator are examined. In particular, a criterion is established for 
determining whether the net effect of seed addition is to augment or decrease the power 
density in MHD generator channels that use magnetically induced nonequilibrium ioniza- 
tion. The analysis considers the effects of electron-neutral collisions, electron-ion 
collisions, ion slip, and variation of the electron energy loss factor with seed fraction 
and seed material. Included in the consideration of ion slip are the effects of seed frac- 
tion and seed material on ion-neutral collision frequency and on the relative drifts be- 
tween seed ions and carrier ions. 

The presented criterion is based on the determination of a break-even condition; 
namely, that combination of gas density, magnetic field, and induced voltage for which 
the power density of the unseeded carrier is neither augmented nor decreased by the 
addition of seed. The formulation is given in general terms and is applicable for arbi- 
trary combinations of monatomic carrier and seed materials. The criterion may be 
separately applied to each successive station of a segmented-electrode Faraday gener- 
ator. 

Calculations of the break-even condition for various noble carrier gases with alkali 
seed have been performed and are presented in graphical form. The results indicate 
that ion slip very strongly affects the break-even conditions. Generally speaking, seed- 
ing is desirable, and the most beneficial effects are obtained for heavy seed 

INTRODUCTION 

One of the key factors in obtaining high power density in magnetohydrodynamic 
(MHD) power generators is the achievement of high electric conductivity in the working 


*The principal results of this investigation were reported at the ASME Winter 
Annual Meeting, New York City, Nov. 29-Dec. 4, 1964. The present report is a con- 
siderably expanded version of material which appeared as ASME Preprint 64-WA/ 
ENER-10 under the title "Concerning the Need for Seeding in Nonequilibrium Magneto- 
hydrodynamic Generator Channels. " 


materials. 




fluid. In closed-cycle electric -power systems where seeded combustion plasmas are 
inapplicable, this problem is particularly severe because, at the highest working fluid 
temperatures currently contemplated for nuclear reactors or heat exchangers, the equi- 
librium electric conductivity of even a seeded gas is far below that required for a rea- 
sonable power density. A nonequilibrium means of obtaining high electric conductivity 
is accordingly of great interest. 

The contemporary studies of nonequilibrium conductivity began with that of Kerre- 
brock (ref. 1). He considered the case of direct current gas discharges and showed that 
the electron temperature, at which the energy received by electrons from the electric 
field is balanced by that lost in collisions, can be considerably above the bulk gas tem- 
perature. He has in fact demonstrated the results to be well represented by a theory 
that assumes the ionization and recombination processes to be in equilibrium at the elec- 
tron temperature so that the electron population is that corresponding to equilibrium at 
the electron temperature. The validity of Kerrebrock's hypothesis has been established 
for a cesium plasma by Ben Daniel and Tamor (ref. 2). 

Hurwitz, Sutton, and Tamor (ref. 3) noted that the electron heating effect might also 
be obtained using the induced electric fields present in MHD generator channels. In 
their paper, as in Kerrebrock's analysis, the electron temperature was determined by 
balancing the energy gained by the electrons from the induced electric field with the en- 
ergy lost by them in elastic collisions with the neutral atoms. They assumed that at the 
low degrees of ionization attained in generators, the electron-ion collisions would occur 
infrequently enough to be considered negligible. 

Reference 4 showed, however, that at the elevated electron temperatures resulting 
from electron heating, a degree of ionization might be attained where the dominant proc- 
ess for momentum and energy exchange is electron-ion collisions. This was found to be 
particularly the case in noble gases displaying the Ramsauer effect. The addition of 
easily ionizable seed to a carrier gas has the two effects of (a) increasing the electron 
density of the working fluid and (b) increasing the electron energy losses through colli- 
sion with seed atoms and ions, thus limiting the electron temperature. The first effect 
tends to augment the electric conductivity, while the second tends to decrease it. It was 
shown in reference 4 that in many cases for argon carrier gas with cesium seed the sec- 
ond effect was predominant, causing an overall reduction in electric conductivity. In 
fact, it was recommended in reference 4 that argon be used without seed as the working 
fluid in an MHD generator. This result is, however, sensitive to the values of the per- 
tinent collision cross sections. The correct calculation of the electrical conductivity 
depends on the average rate of momentum transfer in collisions, which in turn requires 
the averaging of the momentum transfer rate over the velocity distributions of the col- 
liding particles. In the previous investigations (refs. 1, 3, and 4), a value of the mono - 
energetic cross section corresponding to the average electron energy was chosen, which 
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is tantamount to ignoring the distribution of electron velocities. Furthermore, the prior 
investigations did not quantitatively consider the effect of the presence of both seed and 
carrier gas ions on ion slip. 

It is the purpose of the present report to examine in an orderly manner the various 
consequences of the addition of seed to a carrier gas, including the aforementioned 
effects, and to establish a criterion to determine whether seed addition augments or de- 
creases the power density. The formulated criterion may be separately applied to each 
successive station of a segmented electrode Faraday generator. The formulation is 
given in general terms and may be used for arbitrary combinations of monatomic carrier 
and seed materials. 


ANALYSIS 

The criterion mentioned in the INTRODUCTION will now be obtained. The basic 
idea is to determine the condition for which the addition of a vanishingly small amount of 
seed either increases or does not change the power density (power generated per unit 
volume). The condition for which the power density is unc hanged by the addition of seed 
will be called the break-even point. 


Basic Equations 

In reference 4, Ohm's law and the electron energy equation are derived, including 
the effect of ion slip, for a slightly ionized gas containing only one type of ion (ions of the 
carrier gas). In the present report, it is necessary to include both carrier gas ions 
(c + ) and seed ions (s+). (Symbols are defined in appendix A. ) The details of the deriva- 
tion are given in appendix B, where it is shown that all the basic equations of reference 4 
still apply, provided that the definition of the ion-neutral collision frequency v - a is 
appropriately modified (see eq. (B19)). 

The energy balance equation for the electrons is (ref. 4) 


j2 = 3 #( T e - T > N e 


( 1 ) 


where J is the current density, e is the electron charge, k is the Boltzmann constant, 
m is the mass of an atom of the carrier gas, T is the electron temperature, T is 
the static gas temperature, and N g is the number density of electrons. The factor 6, 
which is sometimes called the mean loss factor, is a dimensionless quantity defined for 
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elastic collisions as 
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im y et 
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( 2 ) 


where the subscript t denotes each species of atom or ion that may be present, and v ^ 
is the average collision frequency for momentum transfer between the electrons and 
species t. The collision frequency is written 


y et =N t< v e) Q et 


where ( v ) is the average thermal speed of the electrons 


< v e> 



(3) 


(4) 


and Q et , the average cross section for momentum transfer, is 


Qet< T e> 



v^Q (v )exp 
e mr e' r 


/ m v 2 
- e e 

\ 2kT e 



(t = i, a) 


(Cll) 


The quantity Q m (v e ) appearing in the integrand of equation (Cll) is the momentum 

transfer cross section for a monoenergetic beam of electrons of speed v . Equation 

(Cll) is the result of averaging the rate of momentum transfer in collisions over the 

velocity distributions of both the electrons and the species t of heavy particles (atoms 

or ions) in the manner described in appendix C. When the average momentum transfer 

cross section is calculated from equation (Cll) for any one of the noble gases for some 

ranges of T , the result is found to be as much as an order of magnitude larger than the 

value of CD at beam energies equal to kT Q . The difference between CD (T ) and 

Q is due to the large variation in Q with v for the noble gases displaying the 
m xn c 

Ramsauer effect. (If Q m were independent of v g , this difference would not be signif- 
icant, since in that case Q ec (T e ) = (4/3)Q m . ) Therefore, the use of the monoenergetic 
cross section instead of the average cross section, as was done in references 1, 3, 
and 4, can lead to significant errors. 
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For a segmented-electrode MHD generator, Ohm's law is (ref. 4) 


a (i 

j = H — uB(l - K) (5) 

1 + ^ 

where u is the gas velocity, B is the intensity of the magnetic field (Wb/sq m), and K 
is the load parameter, which is the ratio of the actual voltage across the generator elec- 
trodes to the open-circuit voltage, 



( 6 ) 


The conductivity Oq is 


CT 0 " 


N e e 

m e u e 


(7) 


where 


v e = E v et (8) 

e t/e et 

is the total collision frequency for electrons. The electron Hall parameter /3 g is de- 
fined as 


0 


e 



eB 

"Ve 


(9) 


where o? e is the electron cyclotron frequency. The parameter 0. is 


2eB 


o>i 

^ S ^ =m c v ia. 


( 10 ) 


where v ia is the effective ion-neutral collision frequency. The value of 0 4 is a meas- 
ure of the slip of the ions through the neutral gas. The ion slip varies inversely with the 
gas density (through i> ia ) and directly with the magnetic field. The expression for y- a 
is derived in appendix B (see eq. (B19)). 
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The power density for the segmented electrode generator is 


1 


P = JE = 


N e^ 


3e 2 k6 


m. 


(T - T) KuB 


(ID 


The seed fraction s, defined as 



( 12 ) 


is the parameter of main interest. The number densities of seed and carrier gases in 
the absence of ionization are denoted as N° and N°, respectively, where 

S 


N s = N s + N + 
s 


N? = N c + N 


( 13 ) 

( 14 ) 


The degrees of ionization of the seed and carrier gas are 


respectively. The condition of charge neutrality is 


*This is equivalent to the conventional formula 

a n u 2 B 2 K(l - K) 

P = 

1 + ^ 

Equation (11) is more convenient for subsequent manipulations. 
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isr 


N 




( 15 ) 


( 16 ) 


6 



(17) 


N e * N i - N s+ + N c+ * ( x c + SI X 


It is assumed that the various ion and atom populations are in equilibrium with the 
free electrons at the electron temperature T . Then the Saha equations (eqs. (4) and (5) 
of ref. 4) may be written for the seed as 


( x c + sx s )x s 

1 - x„ 


2g + /2fmkT\ 3 ^ 2 


N° g s 


exp — 


.!L 

kT, 


and for the carrier gas as 


(x„ + sx„)x„ 
v c s' c 


2g + 
_1 (T 

K gc 


/277m e kT e 


3/2 


exp 


kT, 


(18) 


(19) 


where g and g are statistical weights for the ground states of the seed and carrier 
s c 

atoms, respectively, g and g are statistical weights of the ground states of the 

s + c + 

seed and carrier ions, respectively, and V and V_ are the first ionization potentials 

S 

of the seed and carrier gas, respectively. 


Condition for Operation Without Seed 

The power density, current, conductivity, electron density, and temperature all 
depend on the seed fraction s. Suppose that the generator is operating without seed 
with a prescribed load parameter K, magnetic field B, and gas velocity u. The prob- 
lem is to determine the direction in which the power density varies when a small amount 
of seed is added. For this the sign of dP/ds in the neighborhood of s = 0 must be de- 
termined. When dP/ds = 0, the break-even condition will be attained, and this condi- 
tion will require a certain relation between the operating parameters. 

With the aforementioned assumptions, the electric field E is fixed, and the frac- 
tional change in the power density caused by the addition of seed is, from equation (11), 

JL dP = 1 dJ = J_ + 1 1 dT e + J_ d6 ( 2 0) 

P ds J ds N e ds 2 T e - T ds 26 ds 

From equations (1), (5), and (7), it follows that for fixed K, u, and B 
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( 21 ) 


"e (1 + V/ (T e " T)6 = const 


Logarithmic differentiation of equation (21) gives 



1 + Vi 



dy. \ 

1 

ia 1 

y ia 

ds / 


T ds 6 ds 


( 22 ) 


Since v , v., and 6 are known functions of s and T (s), equation (22) may be used to 

G la c 

determine dT 0 /ds. Then the resulting expression for dT /ds can be substituted into 
equation (20) to yield an expression for dP/ds. The details of this derivation are given 
in appendix D, where it is shown that dP/ds i 0 if 



The number n is simply the exponent in a semiempirical relation for the cross section 
Q ec ~ Tg for scattering of electrons by carrier gas atoms. Typically, n is of the 
order 1. The quantity 1 k, which is defined in equation (D18), is related to the differ- 
ence in mobilities of the two ion species in the carrier gas. 

The aforementioned relation (eq. (23)) is derived for the limit of zero seed. It has 
been assumed that as s approaches zero, the degree of ionization of the seed ap- 
proaches one; hence, the ionization potential and cross section of the seed do not appear. 
From equations (18) and (19), it can be shown that for s = 0, x « 1 to a high degree of 
approximation. 

The inequality (23) is an implicit relation between T g , N^, B, and T, which estab- 
lishes a range of electron temperatures in which the power density is not increased by 
seeding. It is not possible to obtain an explicit relation, because the various terms of 
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equation (23) are not simple functions of T . It is convenient for calculational purposes 
to assign N°, T g , and T and to calculate /3 e /3. } which is proportional to B 2 . For 
fixed values of T g and T, the limiting magnetic field is then calculated over the allow- 
able range of gas densities N^. From curves of B against N°, the break-even value 
of corresponding to any given magnetic field may be found for prescribed values of 
T e and T. 

The electron temperature is, of course, not an independent parameter, because it 
depends on the rate of Joule heating through equation (1). Thus it is necessary to estab- 
lish whether or not any prescribed electron temperature is actually attainable. The cri- 
terion for this is the value of the gas velocity u required to produce a given electron 
temperature. The required gas velocity is obtained from the relation 

“BU-K) = ^yf c a + (yy* e> /T e - t <24> 

which follows from equations (1), (5), and (7). (Note that 6=1 in the limit of zero 
seed. ) Finally, the break-even condition is conveniently expressed as a family of 
curves of uB(l - K) against for various values of B and T. In appendix E, a 
sample calculation is given in order to clarify the details of the calculational procedure. 


RESULTS AND DISCUSSION 

Although the break-even condition (the equality in relation (23)) is utilized in the 
manner just described to obtain numerical results, it is of interest to consider the 
simpler forms to which it reduces in special cases so that the physical meaning of the 
condition may become clearer. 


Zero Ion Slip 

A useful limiting case is that in which ion slip is neglected. Then ^ = 0, and the 
right side of relation (23) must also be zero. The break-even condition for zero ion slip 
is further simplified for small degrees of ionization. For this case (x c « 1 and 
s = 0), equation (19) may be written 


g< T e> 

( n £) 1/2 


( 25 ) 


9 


where 


g(T e )- 




From equations (3), (15), and (17), for x £ « 1, 


^ei 


ec 


Q, 


1 - x c %c 


Q • 

ei _ ^ei 
J x c 


Q, 


ec 


(26) 


(27) 


Thus, for « 1 and s = 0, equation (23) yields the following quadratic equation for 
the ratio v ec /v ei : 


a 



( v ec\ 

+ b — - cgO 
Y ei/ 


(28) 


where 


a = 2n + 1 + 



(29a) 


b = 


2n + 1 + 



3 + 





- (2n + 1) 



(29b) 


c 


= 3 



Q, 


ec 


Q, 


ei 



(29c) 


Break-even operation of course corresponds to the equality sign in the relation (28). 
Since a and c are positive, equation (28) has only one positive root, namely, 




+ 4ac - b 
2a 


(30) 


The corresponding gas density is 


10 


( 31 ) 



which is the gas density for break-even operation with zero ion slip. At higher gas 
densities, v ec /^ e i > ^ec^^ei^lim’ inequality (28) is not satisfied. Hence, seeding 
increases the power density. At gas densities below (N^)^ m , seeding will decrease the 
power density. It is not difficult to understand why the ratio of collision frequencies is 
the critical parameter for break-even operation without seed. At gas densities larger 
than that given by equation (31), electron-atom collisions account for most of the re- 
sistivity of the gas. According to equations (3), (7), and (8), the conductivity in that case 
is proportional to the degree of ionization. Seeding would then increase the degree of 
ionization and therefore improve the conductivity and power. At gas densities lower than 
break-even, where i> e - is much larger than u QC , the conductivity becomes independent 
of the electron density and increases with the electron temperature. In this case, the 
addition of seed would lower the electron temperature (cf. eq. (D21)) and consequently 
decrease the conductivity and power density more than could be offset by the ionization 
of the seed. 

Equation (31) sets an upper limit on the gas density, above which the addition of seed 
will increase the power density. Strictly speaking, it is an implicit relation for (N^) 1 im , 
inasmuch as the gas density also enters into Q gi through the term In A (see eqs. 

(C19) and (C20)). Equation (31) is used in the following manner to establish the break- 
even condition for zero ion slip. The electron and gas temperatures are assigned, a 
value for In A between 5 and 10 is assumed, and (N°)^ im is calculated from equa- 
tion (31). The value of In A is recomputed with this density, and the calculation is re- 
peated a few times until it converges. The corresponding value of uB(l - K) is then 
calculated from equation (24). Appendix E contains a sample calculation. 

The results of such a calculation are conveniently represented in a graph of the 
voltage uB(l - K) plotted against the gas density N°. This form is convenient because 
the variation of the gas dynamic properties along a duct can easily be traced. For a 
constant velocity generator for which B and K are also constant, the state of the gas 
would be represented by a horizontal line, while a constant area generator with constant 
B and K would be represented on a log-log plot by a straight line of slope - 1, since pu 
is constant. In general, K is not constant, and the generator operating curve is not 
straight, as will be seen subsequently. 

The zero-ion-slip break-even condition for some of the noble gases is shown in fig- 
ure 1, where the atomic mass of the seed is assumed to be equal to that of the carrier 
gas in order to reduce the number of parameters. For some combinations of gas and 
seed such as argon and potassium or xenon and cesium, this is very nearly the case. 
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Gas pressure, atm 

Figure 1, - Break-even curves for zero ion slip and equal atomic masses of carrier 
gas and seed. Static gas temperature, 500° K. 


When m = m , b = a - c and equation (30) becomes 
s c 



(32) 


- 3 

since Q ec /Q e j is usually of the order 10” . Equation (32) can also be obtained quite 
directly from the right side of equation (23) by setting m = m and x « 0, From 
equation (32), it is clear that the electron-ion and electron-atom collision frequencies 
are of the same order at break-even. 

The break-even curves in figure 1 are for all practical purposes straight lines. 

There is no obvious reason why this should be so, but it is most likely due to the fact 

that the exponential variation of g(T ) in the formula for (N°)„ (eq. (31)) is much more 

e c lim 

significant than the purely algebraic variation of the other terms. To the right of any one 

of these lines, N° > (N°) , and the power density is increased by seeding, whereas 

c c lim 
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Gas pressure, atm 

Figure 2. - Effect of seeding with material whose atomic mass is greater than that of 
carrier gas. Argon; static gas temperature, 5D0° K. 


seeding decreases the power density to the left of the line. 

There is no easily discernible pattern in the location of the break-even curves for 
different gases. It is evident from equations (31) and (26) that for a given electron tem- 
perature (N^)j. will increase exponentially with decreasing ionization potential and 
varies as the inverse square of the cross section Q gc . Furthermore, uB(l - K) is 
proportional to Q , since in the present case v » v, (1 + a/c). Therefore the break- 
even curve for the gas with the lowest ionization potential should be closest to the lower 
right corner of the graph. Unfortunately, for the noble gases a low ionization potential 
does not go hand in hand with a low electron scattering cross section. In fact xenon, 
which has the lowest ionization potential, has one of the highest cross sections. The 
conflicting requirements on the ionization potential and cross section are responsible for 
the location of the curves shown in figure 1. 

The obvious way to use the curves of figure 1 is to choose a particular gas and the 
values of the operating parameters u, B, K, and n£ and to decide whether seeding is 
required by where the point representing the operating condition falls relative to the 
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break-even curve. If the break-even curve is below the operating point, addition of seed 

with atomic mass equal to that of the carrier gas will decrease the power density; 

however, the addition of seed with atomic mass different from that of the carrier gas 

may move the break-even curve to the other side of the operating point, which indicates 

that such seed will increase the power density. Examination of equations (28) to (30) 

shows that if the seed atoms are heavier than the carrier gas atoms, then the ratio 

(v /v ■) and (N°) are decreased and the break-even curve is moved to the left. 

This behavior may be accounted for by the behavior of the elastic loss factor 6 defined 

by equation (2). If m > m , then the addition of seed decreases 6. Figure 2 illus- 

s c 

strates this effect when argon (atomic weight, 39. 944) is seeded with cesium (atomic 
weight, 132.91). 


Effect of Ion Slip 

For /3^ j/ 0 the break-even condition is evaluated numerically, as described pre- 
viously and in appendix E. It is instructive to consider first the simplest case, where 
both the atomic masses and the mobilities of the seed and carrier gas are equal. In this 
case, the loss factor 6 is unchanged by the addition of seed; and the quantity 
representing the difference in ion mobilities, vanishes. The break-even condition 
becomes 



When N° = (N°)_. , as defined in equation (31), the numerator vanishes and 0/3. = 0. 

c c lim 0 

Suppose that for a fixed T . N is continually decreased below (N ) . It is clear 

from equations (25) and (27) that the degree of ionization will then increase, as will the 
ratio v-/v . Therefore, j3 6. will increase with decreasing N°. (The term contain- 
ing x„ in the denominator is generally of less importance for small x . ) Physically, 
the increase in 0 e /3^ is due to an increase in ion slip. Because of the lower gas density, 
ion-neutral collisions become less frequent, and the mechanism for extracting part of 
the mechanical energy from the gas and using it to Joule heat the electrons becomes less 
effective. As shown in figure 3, the effective electric field uB(l - K) required for elec- 
tron heating is everywhere above that for 0 g 0^ = 0, the difference between the curves 
increasing as /3 /3^ increases (see eq. (24)). This difference becomes large enough to 


14 




.001 .01 .1 1 


Gas temperature, atm 

Figure 3. - Effect of ion slip on break-even curve for argon. Static gas temperature, 

500° K-, m s * m c ; equal ion mobilities. 

turn the break-even curve upward at low density. The range of parameters over which 
seeding is advantageous is increased, since dP/ds > 0 everywhere below the break- 
even curve. 

For unequal atomic masses and mobilities of seed and carrier gases, the situation 
changes markedly from that of the preceding example, because the variations of the loss 
factor 5 and of the effective ion-neutral collision frequency v. become appreciable. 

la 

This can be ascertained from the numerical magnitudes of the second and fourth terms 
on the left side of equation (23), since eV /kT » 1. Because the general break-even 
condition is quite sensitive to the ratios of ion mobilities and atomic masses of the seed 
and carrier, as well as to the cross section and ionization potential of the carrier, it is 
difficult to draw general conclusions from equation (23). Each combination of seed and 
carrier must be considered separately. The results are shown in figures 4 to 9 for a 
static gas temperature of 500° K (which is typical of supersonic generators) and for 
various magnetic fields. Figure 4, for cesium-seeded argon, and figure 9, for cesium- 
seeded helium, are of a different character than figures 5 to 8. The break-even curve 
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.001 .01 .1 1 .001 .01 .1 1 .001 .01 .1 1 

Gas pressure, atm 

(a) Magnetic field strength, 1 weber per square (b) Magnetic field strength, 2 webers per square (c) Magnetic field strength, 5 webers per square 
meter. meter. meter. 

Figure 5. - Break-even curves for potassium-seeded argon. Static gas temperature, 500° K. 
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.001 .01 .1 1 .001 .01 .1 1 .001 .01 .1 1 

Gas pressure, atm 

(a) Magnetic field strength, 1 weber per square (b) Magnetic field strength, 2 webers per square (c) Magnetic field strength, 5 webers per square 
meter. meter. meter. 

figure 7. - Break-even curves for sodium-seeded neon. Static gas temperature, 500° K. 
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.001 .01 .1 1 .001 .01 .1 1 

Gas pressure, atm 

(a! Magnetic field strength, 1 weber per square (b) Magnetic field strength, 2 webers per square 
meter. meter. 



of figure 4 has an upper and a lower branch corresponding to high or low electron tem- 
peratures, respectively. In the regions above the upper branch and below the lower 
branch, seeding decreases the power density, while between the branches seeding in- 
creases the power density. The upper branch is similar to the curve of figure 3 and may 
be interpreted physically in the same way. In the region below the lower curve, the 
electron temperature is low and collisions of electrons with carrier gas atoms limit the 
conductivity. The addition of very small amounts of seed changes the conductivity only 
slightly, but because the seed ions are more mobile than ions of the carrier gas, seeding 
causes ion slip to increase more rapidly than the conductivity, which causes the factor 
Uq/( 1 + (3^) in the power density to decrease. Although seeding is not required in the 
region below the lower branch, this region is of little practical significance because of 
the low electron temperature and consequently the small degree of ionization and low 
power density. 

The break-even curves of figures 5 to 8 have also two branches. The upper branch 
corresponds to high electron temperatures and low gas densities (this branch is not 
shown in figs. 5(c), 6(c), and 7, since it falls outside the range of electric fields and 
gas densities of interest). The lower branch, corresponding to higher gas densities and 
lower electron temperatures, is below and to the right of the break-even line for zero 
ion slip. Instead of crossing the zero ion slip line, as in figures 4 and 9, this lower 
branch becomes tangent to it. The different character of the break-even curves of fig- 
ures 5 to 8 from those of figures 4 and 9 is due to the variation in the mobilities of 
carrier and seed ions among the various combinations of carrier and seed. 

Seeding increases the power density in two regions of figures 5 to 8, above the 
upper branch, and to the right of the lower branch. In the first region, the effective 
electric field uB(l - K) is generally too high to be attained in practice. The range of gas 
densities and field strengths in the second region are of more practical interest. 

A typical generator operating curve is indicated by the short dashed curve in fig- 
ures 4(b) and 5(b) for argon. The generator operating conditions are taken from one of 
the cases that was extensively investigated in reference 5. In that reference the one- 
dimensional MHD equations were solved for a constant-area channel, where the nonequi- 
librium conductivities of seeded rare gases were calculated on the same basis as they 
are herein. The electric field E was taken to be constant in reference 5, as required 
by Maxwell’s equation for the assumed one-dimensional situation. Since the gas velocity 
varies along the generator duct, for this case it is convenient to define a new constant 
load parameter as 


K, = — — 
1 u. B 


(34) 
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where u^ is the gas velocity at the entrance section. When this definition is combined 
with conservation of mass, pu = p^u^, the following equation for the operating curve is 
obtained: 


uB(l - K) = u 1 B 




(35) 


The dashed curves shown in figures 4(b) and 5(b) are for an entrance Mach number of 
3. 0 and an entrance stagnation temperature of 4000° R (2222° K), for which the entrance 
gas velocity u 1 is 1316 meters per second. For a magnetic field strength of 2 webers 
per square meter, a stagnation pressure of 1. 64x10 newtons per square meter 
(1. 62 atm) yields the optimum power density (ref. 5). The same conditions are assumed 
herein. In this case it can be seen that the addition of cesium will increase the power 
density except at the very end of the operating line. This conclusion is in agreement 
with reference 5, in which the addition of a very small amount of cesium seed 
(s = 0. 000366) to argon carrier gas was found to yield the optimum total power output. 
Figure 5(b) shows that seeding with potassium is detrimental to the power density 
throughout the generator. Potassium seed was found to be inferior to cesium in refer- 
ence 5, but potassium- seeded argon still produced a higher total power than unseeded 
argon. The mobility of ions, however, was calculated on the basis of the polarizability 
of the carrier gas atoms in reference 5. Such a calculation gives equal mobilities of 
argon and potassium ions in argon, whereas the mobility of argon ions in argon is less 
than that of potassium ions in argon because of charge exchange. Inclusion of charge ex- 
change in this report accounts for the difference in conclusions regarding seeding of 
argon with potassium. 

The upper limit on the gas velocity is 


TABLE I. - MAXIMUM GAS 
VELOCITIES 


Gas 

Stagnation temperature, Tq, °K 

1500 

2000 

2500 

3000 

Maximum gas velocity, m/sec 

Helium 

3947 

4558 

5096 

5582 

Neon 

1758 

2030 

2270 

2486 

Argon 

1250 

1443 

1613 

1767 

Xenon 

689 

796 

890 

975 


u max t /2c p T 0 (36) 

where c n = — — is the constant pressure speci- 
p 2 rn 

c 

fic heat and Tq is the stagnation temperature. 
The values of u„„„ are useful in interpreting 
the results and are given in table I. The maxi- 
mum attainable electric field and voltage may be 
estimated from these values of u m „ v and the 
magnitude of the available magnetic field. At the 
present time, a magnetic field of 2 webers per 
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Gas density, N°, atoms/cu m 

Figure 10. - Effect of gas temperature on break-even curves for cesium-seeded argon. 

Magnetic field strength, 2 webers per square meter. 

square meter (20 000 G) is typical of MHD generators. For this magnetic field and 
K = 0. 5, the maximum attainable field uB(l - K) is equal to u may in magnitude and is 
of the order of 10^ volts per meter. With a very high magnetic field, say B = 10 
webers per square meter (100 000 G), the magnitude of uB(l - K) may reach 10 4 . 

Effect of Static Temperature 

In figures 1 to 9 the gas temperature was taken to be 500° K, a typical value for 
supersonic generator. The effect of changing the gas temperature from 500° to 1500° K 
is shown in figure 10. The break-even curves for zero ion slip are very slightly af- 
fected by the gas temperature. However, the character of the curves for /3. ^ 0 varies 
significantly with gas temperature in the range from 500° to 600° K; and in fact, for gas 
temperatures of 600° K and higher, the break-even curves for cesium-seeded argon have 
the same character as the curves for the other combinations of gas and seed at 500° K 
shown in figures 5 to 8. Since there is insignificant variation of the alkali ion mobilities 
in the noble gases with temperature, the change in the character of the break-even 
curves that accompanies a change in gas temperature is due to the variation of the 
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mobility of the carrier ions with temperature, caused by charge exchange (see the dis- 
cussion in appendix E). 


Effect of Magnetic Field 

Finally, it is interesting to note that the break-even curves for various magnetic 
fields fall rather close together when u(l - K) is plotted against N°/B, as shown in 
figure 11 for a constant gas temperature of 500° K. This fact may be used to estimate 
roughly the location of the break-even curves for various magnetic fields when the shape 
of the curve for one magnetic field is known. 


GENERATOR PERFORMANCE WITH UNSEEDED GASES 

The purpose of the preceding analysis was to establish the regions of operation of a 
generator where the addition of alkali seed to a given noble carrier gas would increase 
the power density. A comparison of various combinations of seed and carrier gases re- 
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Power density, P(1 - K)IK, W/cu m PgPj Electron temperature, 



(c) Power density against effective electric field. 

Figure 12. - Performance parameters as functions of electric 
field for unseeded argon. Magnetic field strength, 2 webers 
per square meter; static gas temperature, 500° K. 



(c) Power density against effective electric field. 

Figure 13. - Performance parameters as functions of 
electric field for unseeded neon. Magnetic field 
strength, 2 webers per square meter; static gas 
temperature, 500° K. 
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Power density, P(1 - K)/K, W/cu m PgPj Electron temperature. 










10 2 10 3 10 4 


Neon 



10 3 10 4 
Helium 

Effective electric field, uBIl - K), V/m 


Figure 16. - Comparison of power densities of unseeded noble gases on basis of equal stag- 
nation temperatures. Magnetic field strength, 2 webers per square meter; static gas 
temperature, 500° K; gas density, 10 24 atoms per cubic meter. 

quires the calculation of the power density for each combination at the operating condi- 
tions of interest and for various seed fractions. Since these calculations have been car- 
ried out in reference 5, including the variation in properties along the generator channel, 
extensive calculations of this nature are unwarranted for the purposes of this report. 
However, the magnitude of the power density for zero seed does provide an idea of the 
best choice of carrier gas. For assigned values of B, N°, T g , and T, uB(l - K) is 
calculated from equation (24), and the result is substituted in equation (11) to yield the 
power density parameter P(1 - K)/K. The electron temperature T , and 

P(1 - K)/K are plotted as a function of the effective electric field uB(l - K) in figures 
12 to 15 for unseeded argon, neon, xenon, and helium. Note that the scale of uB(l - K) 
in the helium curves (fig. 15) extends one decade higher than the scales for the other 
gases. At any specified value of uB(l - K), the power densities for argon equal or ex- 
ceed those of the other gases. The single exception is xenon, which produces higher 
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Figure 17. - Momentum transfer cross sections of noble gases. 


2*3 

power at a gas density of 10 atoms 
per cubic meter because the mobility 
of xenon ions in xenon is smaller 
than that of argon ions in argon. 

A comparison of power densities 
on the basis of equal stagnation tem- 
peratures, rather than equal values 
of uB(l - K), is shown in figure 16. 
The curves shown are simply the 

curves of figures 12(c), 13(c), 14(c), 

24 

and 15(c) for a gas density of 10 
atoms per cubic meter, a magnetic 
field of 2 webers per square meter, 
and a static gas temperature of 
500° K. The scales of uB(l - K) 
are shifted so that a single vertical 
line corresponds to the same stagna- 
tion temperature for all four of the 
carrier gases. The dashed lines in- 
dicate a few stagnation temperatures. 
Any two of the abscissa scales are 
shifted with respect to each other by 
an amount equal to the square root 


of the ratio of their atomic weights, since u = ^ 2c^(Tq - T) and c p °c l/m c . 

Figure 16 discloses that argon produces the highest power at low values of uB(l-K), 
but as uB(l - K) increases, the power production in neon and helium exceeds that of 
argon. Although helium yields the highest power densities when uB(l - K) is greater 
than about 4x10 volts per meter, argon and neon appear to be the best for operation 
over a wide range of uB(l - K). The latter are in fact the carrier gases with the lowest 
electron elastic scattering cross sections at temperatures of the order of 1 electron 
volt. The complete channel calculations of Heighway and Nichols (ref. 5) showed that 
with regard to neon and argon carriers, the optimum overall power output for cesium- 
seeded argon was 15 percent larger than that for cesium- seeded neon. At the entrance 
section of their generator, however, neon-cesium yielded a power density about twice 
that for argon- cesium. The overall power output of cesium- seeded helium was found to 
be considerably less than that of argon and neon. These results can be understood on 
the basis of figure 16, since the power density in neon and helium, although larger at 
high values of uB(l - K), fall off with decreasing uB(l - K) faster than it does in argon. 
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In a supersonic, constant-area generator, uB(l - K) drops off rapidly, as shown by the 
operating curve of figure 4(b), because u decreases and K increases, and therefore 
argon, which produces the highest power density at low values of uB(l - K) , also pro- 
duces the highest overall power output. Although the gas density does not remain strictly 
constant through the generator, its variation is not significant, as can be seen from fig- 
ure 4(b), and the interpretation of the results of reference 5 by means of the constant 
gas density curves of figure 16 should be qualitatively correct. The behavior of the 
curves of figure 16 is not surprising in view of the variation of the elastic- scattering 
cross sections of the noble gases (fig. 17 , p. 26) with decreasing electron temperature. 

The vertical lines in figure 16 corresponding to various stagnation temperatures 
were drawn under the assumption that B(1 - K) = 1 or K = 0. 5 for the particular case 
shown of B = 2 webers per square meter. 


CONCLUDING REMARKS 

The consequences of the addition of alkali seed to a noble carrier gas in an MHD 
generator have been examined from the point of view of realizing magnetically induced 
nonequilibrium ionization. An examination of the physical bases of the results tends to 
identify the desirable properties of carrier and seed materials for this optimization of 
pow 7 er density. 

The important properties of carrier gases are that (1) they be of low molecular 
weight so as to yield large induced electric fields for a given stagnation temperature and 
(2) they have low cross sections to elastic scattering of electrons so that the electrons 
will tend to higher temperatures, resulting in high electric conductivity. The seed on 
the other hand should have (1) low ionization potential so that in low concentrations it is 
fully ionized, (2) high molecular weight to minimize electron energy loss, and (3) low 
cross sections to elastic scattering of electrons for the reason mentioned previously with 
regard to the carrier gas. 

For the gas combinations considered herein, neon and argon appear to be the most 
suitable carrier gases while cesium is the most suitable seed. For reasonable generator 
operating conditions using these carriers, higher power densities will likely be obtained 
with seed rather than in the absence of seed and with cesium seed rather than with any 
lighter alkali gas. 

It should be pointed out that the present analysis directly answers the question of 
whether the power density increases or decreases with the addition of seed only in the 
limiting case where the seed fraction approaches zero. The conclusions cannot be ex- 
tended to arbitrary seed fractions since the power density does not necessarily vary 
monotonically with seed fraction. In fact, the results of reference 5 indicate that for a 
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case where the present analysis indicates the desirability of seeding there is an optimum 
seed fraction beyond which power density decreases with increase in seed fraction. 

Another important restriction of the present analysis is due to the omission of in- 
elastic collisions and radiation losses. The analysis of the combined effect of these 
processes on nonequilibrium ionization in a seeded plasma is a formidable problem, but 
one which is important enough to warrant further study. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 12, 1965. 
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APPENDIX A 


SYMBOLS 


Mks units are used except where explicitly noted otherwise. 


A 

atomic weight 

N 

r 

number of particles or 
th 

a 

B 

defined by eq. (29a) 
magnetic field strength 

K 

r species per unit 
volume 

sum of number densities 

b 

defined by eq. (29b) 


of ions and atoms of 

b o 

impact parameter for 90° 


species r 

Cp C 2 , Cg, C 4 

coulomb scattering 
defined by eqs. (B15) 

N 0 

defined by eq. (31) 
standard gas density, 

c 

defined by eq. (29c) 

IQ 

2. 69x10 atoms/cu cm 

D 

defined by eq. (D22) 

n 

exponent in eq. (Dll) 

d 

E 

Debye length, eq. (C14) 
electric field 

P 

Q 

power density 
total cross section for 

e 

electron charge, ■ 
1.602X10" 19 C 


scattering of mono- 
energetic beam, 

g 

ground state degeneracy 


eq. (C25) 

g(T e ) 

h 

defined by eq. (26) 

Planck constant, 

6. 625x10 3 ^ J-sec 

Qm 

momentum transfer cross 
section for scattering 
of monoenergetic beam 
defined by eq. (C6) 

J 

current density 

Srt 

collision cross section 

K 

load parameter defined 
byeq. (6) 

for momentum transfer 
from species r to 

*rt 

defined by eq. (C3) 


species t averaged 
over Maxwell distribu- 

k 

Boltzmann constant, 
1. 380X10" 23 J/°K 


tion, eq. (CIO) 

m r 

|i. 

mass of particle of r 1 

R . 
rt 

defined by eq. (D13) 


species 

s 

seed fraction defined by 

m rt 

reduced mass of particles 


eq. (12) 


r and t 

T 

static gas temperature 
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th 

T r temperature of r n species 

u gas velocity 

XI- 

V r ionization potential of r Ln 

species 

v relative speed of particles 

th 

v f velocity of particle of r n 

species 

difference in ion drift velocities 
defined by eq. (19) 

mean ion drift velocity defined 
by eq. (B8) 

w r drift velocity of r^ 1 species 

x c degree of ionization of carrier 

gas defined by eq. (15) 

x g degree of ionization of seed de- 
fined by eq. (16) 

Z^ 1 ) cross section for momentum 

transfer from species r to 
species t averaged over 
Maxwell distribution defined 
by eq. (C5) 

a mean thermal speed, eq. (C4) 

|3 electron Hall parameter defined 

by eq. (9) 

/I ion slip parameter defined by 

eq. (10) 

6 mean elastic energy loss factor 

defined by eq. (2) 

e defined by eq. (B31) 

e Q permittivity of free space, 

8. 8 54X10“ 12 F/m 


A ratio of Debye radius to impact 

parameter for 90° scattering 
defined by eq. (C16), eq. (C19) 

X defined by eq. (B20) 

ji ion mobility 

Pq ion mobility at 300° K and gas 

density of 2. 69x10^ 
atoms/cu cm 

v total electron collision frequency, 

eq. (8) 

v- effective ion-atom collision fre- 

la 

quency defined by eq. (B19) 

i; t mean collision frequency for 

momentum transfer from 
species r to species t 

p gas density 

<Tq electrical conductivity 

X scattering angle in center of 

mass coordinates 

defined by eq. (D18) 
n scattering solid angle 

u> e electron cyclotron frequency 

Subscripts: 

a atom 

c carrier gas atom 

c + carrier gas ion 

e electron 

i ion 

max maximum 

min minimum 
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i.L 

r r in species 

(r = e, i, a, c, c + , s, s + ) 

s seed atom 

s + seed ion 


. ,th 

t t species 

(t = e, i, a, c, c + , s, s + ) 

0 stagnation conditions 

1 conditions at generator entrance 
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APPENDIX B 


CALCULATION OF ION SLIP FOR TWO ION SPECIES 


The drift motion of a species r of charged particles in 
tic fields is described by the momentum equation (ref. 4) 


e r (E* + w r x B) = £ ™ rt i, rt (w r " ^ t ) 


(r = e, s + , c + ; 


crossed electric and magne- 
t = e, s + , c + , s, c) (Bl) 


where 


E^E + uXB 


(B2) 


is the effective electric field measured in a coordinate system moving with the mean gas 
velocity u, and m^ is the reduced mass of particles r and t, defined as 


m rt = 


m m, 
r t 


m r + m t 


(B3) 


When written out in full for the electrons e, carrier gas ions c + , and seed ions s + , 
equation (Bl) becomes 


e: 


ec 


w. 


w 


“e* + 
es 


w_ - w 


m e (l, ec 


u es )™ e = " e ( E * + w e x B) 


(B4) 


c+: m eW^c+ - ™e> + K+c^c+c + m c+sW™c+ 

+ m , .v (w - w j = e(E* + w XB) (B5) 

c s c + s \ c s / c 

s+: - 5 e> + KV.tc + n W's+s> 5 s+ 

+ m v /w - w \ = e(E* + w XB) (B6) 

s + c + s c + \ s + c7 s + 
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As in reference 4, the atom drift velocities are neglected, an assumption that is valid 
for small degrees of ionization. 

In reference 4, the drift motion of a single ion species for an unseeded plasma in 
crossed electric and magnetic fields was described by the momentum equation 


m„ 


m e y ie( w i ~ ™e) + ~ y ia*i = e(E* +WjXB) 

a 


(B7) 


The aim here is to derive a similar equation for the combined drift motion of the two ion 
species in the seeded plasma. To accomplish this, it is convenient to define a mean ion 
drift velocity as 


_ N c+^c+ + W 


w i = 


N c + + V 


The difference in ion drift velocities 


w ■ 


w c+ - w s+ 


(B8) 


(B9) 


is also a convenient quantity. 

The following equation is obtained by multiplying equation (B5) by N c+ , equa- 
tion (B6) by N g+ , adding and making use of the definitions of Wj and w d (eqs. (B8) 
and (B9), respectively): 


“e^ie^i - ™e) + 


V 


V 


N c+ + N g+ (m c+c I/ c+c + m c+sW + Np+ + Nr+ (m s+c^s+c 


c+ ' r 


N c+ N g+ 

+ m s + s l/ s+s^ Wi + ^N c+ + N g+ ) 2 [ m c + c V c + c + “c+s^c+s^ 
- K+c^s+c + m s+sW] ^d = e(g * + ^i x5 > ( B1 °) 


In reducing the first term of equation (BIO) to this form, use was made of the fact that 


v , = V A- = V- 

c + e s + e le 


(Bll) 
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since the ion-electron collision frequency involves only the electron number density and 
electron temperature and not the properties of the individual ions (provided all the ions 
are singly charged). Also, ion-ion collisions cancel out because of conservation of 
momentum in collisions. 

An equation for the difference velocity w^ is formed by subtracting equation (B6) 
from equation (B5). The resulting equation becomes 

V 

< m c+s+Vs+ + m s+cVs+c+ + "Vie* + N + N '"VcVc + "Ws+s* 

L N c + + r V 


N 


( m c+c^c+c + m c+s^c+s) 


N+ + N+ ^ c.ccc cs e.sj W d- ew d XB 

= - [< m c+cVc + m cVc+s* - K+c Vc + "Ws+s*] "i < B12 > 


If 


m c+c^c+c + “c+s^c+s = “s+c's+c + “Ws+s 


(B13) 


then the ions drift with the same velocity. Since v ^ is proportional to N t , equa- 
tion (B13) can be satisfied for arbitrary seed and carrier gas densities only if the col- 
lision coefficients for each ion species with the seed and carrier gas atoms are the same, 
that is, only if 


m cVc+c = m sVs+c < Bl4a > 

m cVc+s = m s+sVs < B14b * 


Unfortunately, conditions (B14) are generally not satisfied. Therefore, it is neces- 
sary to solve equation (B12) for w^ in terms of w^ and substitute the result in equa- 
tion (BIO). The following abbreviations are introduced: 

C 1 = m c+c^c+c + m c+s"c+s < B15a ) 



“s+c^s+c + “Ws+s 


(B15b) 
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(B15c) 


^3 m C + S + ^C + S + + m c+r>+ t 'c+f'+ + m o l/ i 


'S T C T S^C -1 


e le 


C 4 = C 3 + 


N c +C 2 + N s+ C l 
N c+ + N s+ 


When compared to the other terms in equation (B15c), m g M ie is of the order (m g 
and therefore may be neglected. Then equation (B12) may be rewritten as 



which may be easily solved by premultiplying both sides by the operator 



and using the fact that is perpendicular to B. Then 



If equation (B17) is substituted into equation (BIO), the resulting equation can be 
written in a form quite similar to that of equation (B7), namely, 


m - - 

m e v ie&i v ia™i = e [5* + (i + Ajw* x 5] 


where now the effective ion-neutral collision frequency is defined as 


^ia 


2 

x c 

m„ 

X„ 

+ sx„ 

c 

c 

s 


sx_ 


C l + 


sx s x c < C 1 - C 2^ 2 1 


x c + sx s (x c + sx g ) 2 C 4 1 + ( — \ 2 


and 


(B15d) 

/m A ) 

(B16) 

(B17) 

(B18) 

(B19) 
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sx s x c ( C 1 ~ C 2) 2 1_ 

(x c + sx g ) 2 cl 1+ /eB 


(B20) 


The collision frequency for electrons with any single ion species A + in a mixture of 
singly charged ions is 

■W = N A+ Q ei< v e> (B21 

where the electron-ion collision cross section is determined by the coulomb inter- 
action between electron and ion and is therefore the same for all singly charged ions 
(see appendix C for the exact expression). In the present case, the total electron-ion 
collision frequency is 

v el = v ec+ + "es+ ' N i Q ei< v e> <B22 


where 


N i " N c+ + N s+ 


If the total electron- atom collision frequency is written 


v = v + p 
ea ec es 


and use is made of the definition of w^ (eq. (B8)) and equations (B22) and (B23), equa- 
tion (B4) may be written in the form 

m e v ei&e ' + "Vea^e = ~ e( ^* + X < B2E 

which is identical to equation (36) of reference 4. 

Ohm’s law and the expression for the electron heating rate are obtained from equa- 
tions (B18) and (B25). The total current 

J = - Wg) (B2( 


and the electron current 
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J e * - N e ew e 


(B27) 


are introduced into equations (B25) and (B18), which become 

-m e ^ el J - m e „ ea J e - -N e e 2 E* + eJ e x B 


(B28) 


— m 


m 


e^ei^ = N p e 2 I* + e(l + A)(J -J p )XB 


la 


(B29) 


Addition of equations (B28) and (B29) and solution of the resulting equation for J g by 
the same device as employed in the solution of equation (B16) yield 


J = 
e 


A — {[l + e + Ml 1- A)d?]j + (ill + e(l + A)1 5 X j] 

(1 + ef +A t J B J 


(B30) 


where 


e = 2 


in 

e ea 

in v- 
c xa 


(B31) 


and 


*i a 


2eB 


in ia 
c la 


( 10 ) 


Equation (B30) is substituted into equation (B28), which may then be put in the following 
form: 


o E* = (i^l T(i + e ) 2 + * 2 e 2 l 

° <l + e) 2 + A 2 fi?\U L lJ "e L 


1 + c + A(1 + A)A 


6 


+ /3 e /3. [l + e (1 + A)]| J - /3 g jl + e + A(1 + A)0 2 - e |l + e (1 + A)j|| Xj) (B32) 
where Oq is the conductivity 
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(?) 


CT 0 = 


Nel 

m t* 
e e 


and the electron Hall parameter /3 p is defined as 

eB 




^e = 


v e ‘Ve 


( 9 ) 


The electron heating rate J e * E* is, from equations (B30) and (B32) 


— — 

J e • E*=i- 

^ rr 


[< 


1 +e) 2 + A 2 ^ 2 ] 2 ^ 


1 + e + A(1 + A)/3 f j(l + e ) 2 + A 2 /3f 


f] 


+ |^1 + 6 + A(1 + + PqP i^ £l + e (1 + A)j y (B33) 


0 

Since e and /3. are both much smaller than 1, and A for small seed fractions 

1 2 2 

(s < 0. 1) is of order 1 or less, terms of order e, eA, and A j3. can be neglected in 

equations (B32) and (B33). Thus, the simplified equations are 






(B34) 



(B35) 


which are formally identical to equations (49) and (51) of reference 4. The only modifi- 
cation is that v- is defined differently now (eq. (B19)). In the limit of zero seed, v. 

ia la 

reduces to v Q+c , as it should. 
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APPENDIX C 


MOMENTUM TRANSFER CROSS SECTIONS AND COLLISION FREQUENCIES 

The collision frequencies ^ used in this report refer to the average rate of 
momentum transfer from species r to species t per particle of species r. The 
average rate of momentum transfer from species r is 



(Cl) 


where ( 6 f r / 6 t) co ^ denotes the rate of change of the distribution function f r due to col- 
lisions. Burgers (refs. 6 and 7) has obtained an explicit expression for (Cl) by assum- 
ing that the distribution function is close to Maxwellian. In the first approximation, he 
obtains 



(C2) 


where w and w, are the drift velocities. The friction coefficient K . is defined as 
r t rt 


^rt 3 m rt N r N t aZ rt 


( 11 ) 


where m^ is the reduced mass 


m rt ~ 


m m, 
r t 


m r + m t 


(C3) 


(B3) 


a 


is a mean thermal speed 


2 


^ As noted in reference 4, Burgers* results apply to only the case T r = T t = T, but 
equation (C4) represents the proper generalization of Burgers’ work to the case of dif- 
ferent temperatures. 
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where v = |v - v.| is the relative speed. 

The quantity Q m (which Burgers denotes S^) appearing in the integrand is the 
momentum transfer collision cross section for a beam of particles of velocity v and is 
defined as 


Q m (v) s - cos y)df2 


(C6) 


where n(v, y) is the differential cross section for scattering through the angle y and 
dfl the element of solid angle. 

The relation between v ^ and K rt is (refs. 4 and 8) 


K 


V-i = 


rt 


rt "Wr 


= - QfN.Z 
3 1 


( 11 ) 

rt 


(C7) 


More often, the 2/3 factor is absorbed into the cross section, and the average value of 
the relative speed, which is 


/8kT 8k T. 

< v ) =[ - r + 

7Tm r 7rm^ 


, 1/2 


(C8) 


is used in place of a. Then is written in the conventional definition as 


v rt = N t< v )Qrt 


(C9) 


Clearly, 
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(CIO) 



( 11 ) 

rt 


For collisions between electrons and atoms or ions, the relative speed v is very 
nearly equal to v , and a « (2kT e /m e )^^; hence equations (CIO) and (C5) may be re- 
duced to 


Q et (T e) 3 1 2kT 



V nQ m ( V J eX P I- 


e^m' e 


2 ' 
m v„ 
e e 

2kT 


dv 


(t = i, a) 


(Cll) 


The well-known Rutherford cross section 


CT (v e , x) = 


/ 2 
4 \4ireJ 


1 

2 4 


m v„ 
e e 



(C12) 


is used to calculate Q m for collisions of electrons with singly charged ions . When 
equation (C12) and the formula df2 = 2n sin x dx are substituted into equation (C6) and 
the integration is performed over all scattering angles from 0 to it, the result is 



(CIS) 


Due to the long range nature of the coulomb potential, this integral diverges at the lower 
limit. Actually, the max imum range of the electrostatic field of an ion is of the order of 
the Debye length, which is (ref. 9) 



(C14) 


At distances larger than d, electrons almost completely screen the field of the ion. 
Hence it is customary to replace the lower limit in equation (C13) by the minimum scat- 
tering angle x m j n ) which is determined from the relation 
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tan 


^min _ _1 
2 _ A 


(C15) 


where A is the ratio of the Debye length to the impact parameter b Q for 90 scattering 
(ref. 9) 


V T e 


, 1/2 


a4 = 

b o 


l Ng e 


2 ! 
e 1 


(C16) 


477e o m v 2 
e e 


Ions also contribute to the screening effect, but as pointed out in reference 10, the time 
involved in an encounter between an electron and an ion is much smaller than the period 
of oscillation of the ions; hence, there is not enough time for the ions to respond to the 
fluctuation in charge caused by the passing electron, and the ions therefore cannot pro- 
vide effective screening for electron-ion collisions. 

The result of cutting off the integral in Q m at x m j n I s therefore 



2t r 



1 

2 4 


in v 
e e 


ln(A 2 + 1) 


or since 


A 2 » 1 


(C17) 



4 3T 



1 

2 4 


rn v 
e e 


In A 


(CIS) 


On account of the slow variation of In A with v , it is also customary to replace ni e v e 
in equation (C16) by its average value, 3kT g , so that A becomes 


A = 


l 2 ^o kT e> 

eV / 2 


3/2 


(C19) 
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Integration of the expression (C18) according to equation (Cll) yields finally the average 
momentum transfer cross section for electron-ion collisions 




In A 
(kT e ) 2 


(C20) 


The ion-ion cross section for singly charged ions is derived from equations (C5) 
and (CIO) by using a Rutherford cross section similar to equation (C12), but with m g 
replaced by the reduced mass of the ions and v g by the relative speed v of the ions. 

The screening due to ions must also be included, since the ion-ion collision time is of 
the same order as the ion oscillation period. The correct Debye length d f , including the 
effect of ion screening, is (ref. 10) 



where 




(C21) 


(C22) 


The result for the ion-ion cross section is 



2 

£ / e 2 \ ln(A ,2 + 1) 
3 \ 47re o/ (kTj) 2 


(C23) 


where A’ is the ratio of the Debye length d’ to the average impact parameter for 90° 

2 

scattering of ions (e /47re 0 )/3kT| 


12*(e 0 kT.) 3/ * 

eV / 2 



(C24) 


The ion temperature is assumed to be equal to the gas temperature T. Since 
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T g » T, the last factor in equation (C24) may be set equal to one. However, A’ « A, 
and it is generally incorrect to make the approximation A* » 1 in equation (C23). 

For electron-atom collisions, experimental values of the momentum transfer cross 
section Q m were used. The authors of reference 5 have evaluated Q e j-(T e ) numerically 
for a number of gases in the range of T g from 100° to 15 000° K. Where experimental 
values of Q m were not available, they used the total cross section 

Q c (v e )= / CT ( v e’X)<^ ( c25 ) 

instead. The difference between Q m and Q c for most gases is less than 10 percent 
(ref. 11, p. 31). The results for Q ec (T e ) for the noble gases are summarized in fig- 
ure 17 on page 26. 

In reference 4, it was assumed that electron-atom collisions could be treated by the 
hard-sphere model, in which Q m is independent of v . In that case, equation (Cll) 
yields 

Qec = | Qm < c26 > 

- 1 7 

A constant value of Q m = 2. 0x10 square centimeter was used for argon in refer- 
ences 1, 3, and 4. This value was inferred from monoenergetic cross section data 
(ref. 11, p. 19), which showed little variation over a range of electron energies below 
1 electron volt. From figure 17 it can be seen that this value was far too low. The 
average momentum transfer cross section also rises rather rapidly with electron tem- 
perature. The necessity of using the properly averaged cross section is apparent. 

The momentum transfer frequencies for ion-atom collisions are calculated from 
ionic mobility data. The mobility of an ion of species r in a gas of species t is 
defined as 


w r = M rt E (C27) 

where is the drift velocity and E is the electric field. Comparison of equations 
(Bl) and (C27) shows that the mobility for singly charged ions is related to the average 
momentum transfer v ^ by the formula 


M rt = 


Vrt 


(C28) 
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Note that v ^ is proportional and p rt is inversely proportional to the density N^. Ex- 
perimental mobilities are often referred to the standard gas density Nq = 2. 69x10*9 
atoms per cubic centimeter and the corresponding mobility denoted by Pq, as in refer- 
ence 11, so that the mobility p at any other gas density N is 



APPENDIX D 


DERIVATION OF BREAK-EVEN CONDITION 


The derivatives dN /ds, dv /ds, dy. /ds, and d5/ds are required for the evalua- 

tion of dT e /ds and dP/ds from equations (20) and (22). Since the former quantities 

are functions of the degrees of ionization x and x of the carrier and seed gases, it is 

convenient to first calculate dx /ds and dx /ds. 

c s 

Differentiation of equations (18) and (19) with respect to s and solution of the re- 
sulting equations for dx /ds and dx /ds yield 

c s 


dx. 


X c (l - x c ) 


ds x c (l - x c ) + sx g (l - x g ) + x c + sx g 


, .. 3 

(X C + SXs ^2 + kT 


+ sx (1 - xj 


e < V c - V 


S' S' k T 


1 dT e 


T ds s 
e 


■x (Dl) 


r- 


dx 


x s (1 - x s) 


ds x„(l - x„) + sx fl - x„) + x + sx„ 
c v c' s v s' c s 


<x c s’lj ^T~ 


Vi- 


x f .(l - xj 


6 <V C - V s ) 


c' c 7 kT 


1 dT e 

X„ 

T ds s 

J 


(D2) 


Only the limiting forms of these expressions for s = 0 will be used here. It will be as- 
sumed that, in the limit of zero seed, the degree of ionization x of the seed approaches 
unity. Numerical calculations carried out in the course of the work reported in refer- 
ence 4 for s < 0. 01 support this assumption over the range of electron temperatures of 
interest. In fact, if s is set equal to zero in equations (18) and (19) for x £ « 1, the 

following expressions are obtained for x and x : 

c s 

X = Vf~ ( D3a ) 

c v c 


V?e 


(D3b) 
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where F amd F are the right sides of equations (18) and (19), respectively. The 

ionization potential of the carrier gas is several times larger than that of the seed, so 

that F_ « F ; and since F > 1 for the range of electron temperatures of interest 
c s s 

here, x = 1 to a high degree of approximation, 
s 

With the aforementioned assumptions, in the limit of zero seed, equation (Dl) re- 
duces to 



1 - 


2 - 



(D4) 


Since dx /ds appears everywhere multiplied by s, its limiting expression for s = 0 is 

O 

not required. 

From equation (17), for x = 1, 

s 





i_ 

ds Js=0 M 1 - x c>. 


CD5) 


Substitution of equation (D5) into equation (20) yields the expression 


J, /dP\ _1 

P\^Lo =2 


2(1 - xj U eV 


V (2 + _ £ \ +- -e 




2 ” x c \ 2 kT e/ T e- T 


T e\ ds 


1 / d6 > 
+ — 


s=0 x c( 2 - x c> 6 \ ds /s=0| 


(D6) 


The various derivatives in equation (22) will now be calculated. The total electron 
collision frequency is written 


V — V + V + V . 

e ec es ei 


(D7) 


where 


(D8) 

(D9) 


^ec = N c Q ec< v e> = N ? (1 ' x cKc< v e> 
v es = N sQes< v e> = N c s(1 " x s> Q es< v e> 
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v ■ = v + + v + = N.Q . / v \ = N°te + sx )Q . ( v ) 
ei ec + es + i*ei' e' c' c s'^ei' e' 


(DIO) 


_2 

The electron-ion cross section Q gi is proportional to T g , as can be seen from equa- 
tion (C20). For the electron-atom cross sections the variation with electron temperature 
is expressed by relations of the form 


Q ec = (const)T^ 


(Dll) 


in the immediate vicinity of the electron temperature corresponding to zero seed. The 
exponent n is taken from the slope of a log-log plot of Q ec as a function of T g . The 
method of deducing the average momentum transfer cross section Q ec (T e ) from experi- 
mental data on the momentum transfer cross section Q m ( v e ) f° r a monoenergetic beam 
of electrons is described in appendix C. 

The derivative of v with respect to s can be calculated in a straightforward man- 
ner using equations (D7) to (Dll) and (D4). The slight variation of In A with s can be 
neglected. In the limit of zero seed the result becomes 


1 



s=0 



2x c 
2 ' x c 


! + ^j| 

2 kT e/ 



-3 


2 (1 -x c) 

, ev 4 

1 ( 

MT e) 

+ 1 

»-*c ' 

'tO | 
H 

** 

.fD 

T e 

V ds ) 

2 - x 
s=0 x c 



(D12) 


The calculation of dv. /ds is also straightforward. It is convenient to use the ab- 

13 . 

breviation 


m rt^rt 

R r,= -V^ = R tr < Dl3 > 

1N t 

where ion-ion and ion-atom collisions are involved, since this quantity is independent of 
s. Then the quantities C^, Cg, and Cg, which enter into the expression (B19) for ia & 
become 

C 1 = [(! - x c> R c+c + s(1 " x s) R c+s] N c (D14) 
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C 2 = [d - *c>Vc + =<‘ - 

c 3 = < x c + ^rtnS? 


(D15) 

(D16) 


Substitution of equations (D14) to (D16) into equation (B19) and differentiation of the re- 
sult with respect to s yields in the limit of zero seed 


1 1 

/diA\ 

I ia\ 


Ids J 


s=0 



* = - 


s=0 



(D17) 


where the expression (D4) has been used for (dx c /ds) , and the following notation has 
been introduced for brevity: 


R s + c _ A ~ x c^ R c + c + x c R s + c + 

R c+c /(‘-WWc* 


(D18) 


In obtaining equations (D17) and (D18), use was made of the approximation 
2 

(eB/C 4 ) « 1. This approximation is consistent with the assumptions of appendix B 

“ n ft 

since the quantity (eB/C^) is of the same order as Pf. 

Finally, the derivative d6/ds is calculated from the expression 
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1 


in 


'ec + v ec + + m ( 
s 


^es + v es + ) 



(D 19) 


and the result in the limit s = 0 is 


1 

6 


d6\ 

ds /s=0 



(D20) 


Substitution of equations (D12), 
ing equation for (dT ( /ds) g _Q: 


(D17), and (D20) into equation (22) yields the follow- 
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(D21) 
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Substituting equation (D21) into equation (D6) and rearranging yields the result 
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where the numerator N is defined as 
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(D22) 


(D23) 



(D24) 
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For the eases of interest in this report, the denominator D always turns out to be 
positive. Thus the sign of dP/ds is governed by the sign of N. For N §0, seeding 
does not improve the power density. 
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APPENDIX E 


DETAILS OF THE CALCULATION 


The break-even condition (eq. (23)) is rewritten as 


0 0 . 

e i it - T 




2eV T 
(1 - x_)(3 + ^ + 
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3 eV c\^ei 
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2 kT„ \ 2 1 


(El) 


where the magnetic field B enters only through 0 g 0j, which for s = 0 becomes 

eB \ / 2eB 


PJ, = 


e i 


im eV m cVci 


(E2) 


The expressions for ^ ec /^ e and v ^/v in the limit s = 0 are obtained from equa- 
tions (D8) and (DIO) 


v (1 - x )Q 

ec ' c /V£ ec 


v e ^ " x c^ec + x c^ei 


^ei 


x Q . 
c^ei 
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where from equations (C19) and (C20) 


(E3) 
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Q , = 4. 34X10" 18 - ln A (sq m) 
el / \2 
/ kT 


(E5) 


A = 1. 55x10 


13 1 


i^kT X 3 / 2 


N i 


1/2 1 


(E6) 


and the limit of zero seed 


N. = xN° 
1 c c 


(E7) 


The degree of ionization x c of the carrier gas is calculated from equation (19), 
which for s = 0 becomes 


6. 04X10 27 / g c+V kT e 


3/2 


1 - x„ 
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(E8) 


When x « 1, a good approximation to x„ is 

v C 


x c = 


g(T e ) 


7. 77X10 1,5 / g c+ 
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(E9) 


Finally, the ingredients will be assembled for the quantity which was defined in 
equation (D18). The quantities R g+c and R c + c for ion-atom collisions are calculated 
from data on ion mobilities. From the definition (D13) of R^ and equations (C28) and 
(C29), it follows that 


Urt 


“rt N t (^o) rt N 0 


(E10) 


25 

where Nq = 2. 69x10 atoms per cubic meter is the standard gas density and (fig) 
the corresponding mobility. If fig is given in square centimeters per volt per second, 
as is usually the case, 
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R rt = 5. 96X10 


-41 1 


(Ma) 


(Ell) 


0 ; rt 


For collisions between ions of species A + and B + , equations (C9) and (CIO) give 


V A + B + ~ N B +Q ii 


8 kT 

7 rm A + B + 


1/2 


(E12) 


if it is assumed that both ion species are at the temperature T of the gas. When equa- 
tion (E12) and the expression for (eq. (C23)) are substituted into the definition for 
the ion-ion collision coefficient R^+b + (eq. (D13)), the following expression is obtained: 


td i , 0vin -27 1/2 ln(A z + 1) 

R A + B + " 1- 38X10 m A /+ B + 3/2 
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Then 4* may be written 
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and A and A are the atomic weights of carrier gas and seed, respectively. The 

V S 

units of (p 0 ) c+c are square centimeters per volt per second in equation (E16). Note 


54 



TABLE H. - PROPERTIES OF NOBLE GASES 


Gas 

Atomic 
weight of 
carrier 
gas, 

A c 

Ionization 

potential, 

V c> 

eV 

Ground state 
statistical 
weight of 
carrier gas - 

Ion mobility at 300° K and 
IQ 

gas density of 2. 69xlO A 
atoms per cubic centi- 
meter (refs. 12 and 14), 
Mqj 

sq cm/(V)(sec) 

Atom, 

Sc 

Ion, 

V 

Helium 

4.003 

24. 580 

1 

2 

10.8 

Neon 

20. 183 

21.559 

1 

6 

4.2 

Argon 

39. 944 

15. 755 

1 

6 

1.60 

Krypton 

83.80 

13. 996 

1 

6 

.90 

Xenon 

131.30 

12. 127 

1 

6 

.58 


TABLE m. - ATOMIC WEIGHTS AND MOBILITIES 
OF ALKALI SEED IONS IN NOBLE GASES 
[Mobility data taken from table 9-2-2 of ref. 12.] 


Seed 

Atomic 

weight 

of 

seed, 

A s 

Background gas 

Helium 

Neon 

Argon 

Krypton 

Xenon 

Ion mobility at 290° K and gas density of 
1 Q 

2. 69X10 atoms per cubic centimeter, 
^0’ 

sq cm/(V)(sec) 

Lithium 

6.94 

24.4 

11. 1 

4.7 

3.7 

2.8 

Sodium 

22. 991 

22.8 

8.2 

3.0 

2.2 

1.7 

Potassium 

39. 100 

21.6 

7.6 

2.6 

1.9 

1.4 

Rubidium 

85.48 

20.2 

6.8 

2.2 

1.5 

1.0 

Cesium 

132.91 

18.5 

6. 1 

2.1 

1.3 

.91 


that vanishes when the mobilities of the seed and carrier gas ions are equal. 

With the previous expressions it is a simple matter to calculate /3 e /3. from equa- 
tion (El) for prescribed values of T , T, and N°. Then the corresponding value of B 
is calculated from equation (E2), which upon insertion of the expression for v c + c in 
terms of P c + c (eq. (C28)) and the relation between M c + c and the standard mobility 
(M 0 ) c+c (eq. (C29)) becomes 


2 oo U - XcK^e^i 

i 2 = 2. 12X10' 33 c e e 1 

^O^c+c 


(E17) 
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(E18) 


v 
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6. 69X10' 



(1 - x )Q + x Q . 
v c'^ec c^ei 


‘] N ° 


The gas velocity required to attain the prescribed electron temperature is calculated 
from equation (24), which may be written 

uB(l - K) = - 8 -- 98 * 1 ° 10 (1 + ^i)"e (T e " T ) l/2 (E19 > 

a 

A c 

The pertinent properties of the rare gases are summarized in table II. Table El 
gives the atomic weights and mobilities of the alkali seed ions in the noble gases. The 
mobilities of the alkali ions in the noble gases are well predicted by classical collision 
theory in which the interaction of the ion and atom is assumed to be a polarization attrac- 
tion (ref. 12). In this case the theoretical mobility is independent of gas temperature, as 
will be assumed here. The available experimental data (ref. 13) indicate that this is a 
reasonably good assumption. The mobility of the noble gas ions in their parent gases, 
however, decreases slowly with increasing temperature because of the effect of charge 
exchange (ref. 14). The following theoretical expressions for the temperature depen- 
dence are given in reference 14 and are in good agreement with experiment: 

— (He+ in He) = 2. 96X10" 3 T 1 / 2 + 3. llXlO" 2 + 2. 11X10" 2 T' 1 + 6. 66 T' 2 (E20a) 
M 0 

— (Ne+ in Ne) = 8. 69X10" 3 T 1 / 2 + 9. 16X10" 2 + 0. 20 T" 1 + 59. 7 T" 2 (E20b) 

M 0 

— (Ar+ in Ar) = 2. 08x10" 2 T 1 / 2 + 0. 24 + 1. 33 T" 1 + 1. 43xl0 2 T“ 2 (E20c) 

M 0 

These expressions, which are claimed to be valid for 50° to 1200° K, are used in the 
calculations. Unfortunately, no expression for the temperature dependence of xenon was 
given. For xenon the same temperature dependence as for argon is assumed (eq. 

(E20c)), but with the constant 0. 24 replaced by 1. 36 in order to give the correct value of 
M 0 at 300° K. 

To illustrate the procedure, a sample calculation for argon seeded with cesium 
(atomic weight = 132.91) is given. An electron temperature T g = 8000° K (0. 689 ev) 
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and a gas temperature of T = 500° K are chosen. First the break-even point for zero 
ion slip is determined from equation (31). From equation (E9) 


g(T e ) = 7. 77X10 13 V6 (0. 689) 3 / 4 exp 


15. 755 
2(0. 689) 


1. 57X10 9 


Let In A = 7. 0 as a first approximation. Then from equation (E5) 


_ _ (4, 34x10' 18 )7.0 

^ei o 

(0. 689; 


6. 40X10 


-17 


(sq m) 


-20 

The electron-atom cross section for argon at this temperature is Q = 4. 25x10 

o GC 

square meter, so that Q ec /Q e j = 0. 664x10” . The coefficients of equation (28) are 

a = 4. 33, b = 16. 06, c = 0. 935 (for n = 1. 13), and from equation (30) 

= 0. 0582. The limiting gas density is from equation (31) 
ec ei lim 2? _ ^ 

(NPiim = 1. 895x10 m , and for this gas density equation (E6) gives In A = 6. 40. The 

calculation is repeated with this new value of In A. The ratio (v / 1/ .),. is essen- 

n 22 ^ lim 

tially unchanged, while the new value of N5 is 1. 582x10 , for which In A = 6. 44. 

Another iteration gives = 1- 60x10^ and In A = 6. 45, and further iterations are 

not required. 

The value of (N°), . determines the point where P J3. = 0. For other gas den- 
sities, P 8. is calculated from equation (El). For N“ < (Nz),. , both sides of equa- 
tion (El) are positive. For the right side of equation (El) is negative, 

but as the gas density is increased, the coefficient of the left side changes from positive 
to negative. In the present example, this occurs at = 4. 69x10 per cubic meter. 
Since P P. must be positive, for the range of densities 1. 6xl0 22 < N° < 4. 69X10 22 the 
break-even condition cannot be satisfied; hence seed increases the power. For 

> 4. 69X10 22 , the coefficient on the left side of equation (El) is negative, and a posi- 
tive value of /3 g |3. may be found that satisfies the equality. Figure 18(a) shows the 
variation of P with N°. 

The sign of dP/ds in regions off the break-even curve is determined from equa- 
tion (23) as follows. If both the right and left sides of equation (23) are positive and if 
P e @i i s smaller than its break-even value (determined from eq. (El)), then inequality 
(23) is satisfied and dP/ds <0. If P Q P^ is larger than its break-even value, the in- 
equality cannot be satisfied and dP/ds > 0. Therefore, in figure 18(a) dP/ds is positive 
above the left branch of the break-even curve and negative below the curve. When both 
the left and right sides of equation (23) are negative, as in the present example for 

N° > 4. 69X10 22 then for P P- < (P P-) the inequality is not satisfied and 

c ’ e 1 ' e l'break-even 
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(a) PgPj against gas density. (b) Magnetic field against gas density. (c) Effective electric field against gas 

density. 


Figure 18. - Results of sample calculation for cesium-seeded argon. Electron temperature, 8000° K; gas temperature, 500° K. 

dP/ds > O. In the region 1. 6x10^ < N° < 4. 69 x 10 ^, dP/ds > 0, as mentioned pre- 
viously. 

Once /3 g /T has been calculated for assigned values of N°, T , and T, the cor- 
responding values of the magnetic field B and the effective electric field uB(l - K) are 
obtained from equations (E17) and (E19), respectively. The results for the example are 
shown in figures 18(b) and (c). From figure 18(b) the values of N° corresponding to a 
particular magnetic field may be picked off; for example, for B = 2 webers per square 

meter, N° is 5. 8x10^ and 1. 15x10^ per cubic meter. The values of uB(l - K) at 

^ 2 2 
these densities are read from figure 18(c) and are 7. 3x10 and 5. 25x10 volts per meter, 

respectively. These two points lie on the lower branch of the break-even curve of fig- 
ure 4(b). Curves similar to those of figures 18(b) and (c) are constructed for each as- 
signed electron temperature, and the densities N° and the electric field uB(l - K) are 
read from the curves for a specified magnetic field. The results for cesium-seeded 
argon are plotted in figure 4 for various magnetic fields. At high electron temperatures 
(above 9000° K), points on the upper branches of the break-even curves of figure 4 are 
obtained. 
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